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1 Introduction 


HISTORY: The subje ct of this grant, which was originally concerned with 
the study of semiconductor alloys and heterost incut ures of use for pho- 
tovoltaic applications was changed in the last year (HMH) 1o a study 
of metal/graphite bonding in view of the interest in this subject in the 
Materials Branch (Dr. Stephen Pepper) at NASA Lewis Research ( en- 
ter. Progress reports on the previous year's activities were given in the 
renewal proposals. 

MOTIVATION: The study was focused on ( 'u/graphite bonding in view ot 
its practical importance in graphite fiber enforcement of ( 'u lor applica- 
tions where high t hermal conductivity is important. The technological 
problem is that the interfacial bonding between Cu and graphite libers 
is extremely weak. 

PREVIOUS WORK: Wetting studies of the graphite basal plane at 
NASA LeRC by De Vincent et al. [1] by Cu and Cu with various 
alloy additions confirmed this and also indicated that small amounts 
of Ti and Cr could lead to better bonding. 

GOALS: The goals of this proposal were: 

1 . To provide a fundamental science basis lor why 1 he bonding ol ( 'u 
to graphtie is weak. 

2 . To critically evaluate the previous analysis ot the wetting studies, 
in particular with regards to the values used for surface enrgies of 
( 'u and graphite. 

To make recommendations for future experiments or other stud- 
ies which could advance the understanding and solution ol this 
technological problem. 

2 Approach 

hirst-principles electronic structure calculations were used to study ( lie prob- 
lem. These are based on density fucntional theory in the local density ap- 
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proximal ion [2] and l lie use of the linear mulfin-tin orbital band-struct ure 
met hod [2] . 

Calculations were performed lor graphite monolayers, single crystal 
graphite' with the hexagonal AB stacking, bulk Cm (hi {111} surface and 
( hi/graphite superlattices. The st udy is limited to the basal plain' ol graphite 
because this is the graphite plane exposed to Cu in the graphite libers. 

A crilieal study was made of the available' literature on Cu and graphite 
surface energies and combined with the measured contact angles to evaluate 
the experimental adhesion energy. 


3 Technical details and results 

For t his initial study, we attempted to use t he linear muflin-t in orbital met hod 
(LMTO) in the so-called atomic sphere approximation (ASA). I his is the 
simplest version of the method which has the advantage ol being computa- 
tionally cost; effective and easy to use. The disadvantage is that the approx- 
imations made limit the accuracy for certain types ol problems. Unfortu- 
nately. severe limitations due to the method were discovered for t lie present 
problem. Alt hough qualitat ice conclusions about the electronic si met lire and 
insights in the origin of tin 3 weak bonding could be obtained, it was impos- 
sible to obtain quantitative results for adhesion or surface energies lor the 
system under sty tidy using this com pul at ional approach. This problem could 
he solved by using the so-called full-potential LMTO method. ] he time and 
resource's available for this grant, however, did not- allow us to carry out- 
1 liese more demanding calculat ions. Instead, we decided to provide some 
semi-quant it ative results using literature data on suriace energies and the 
measured contact angles of the wetting studies. We also discuss some ol 
the literature available on metal/graphite interaction. Below, we first con- 
sider the results of the ASA-LMTO calculations and next, the results of our 
literat ure studies. 

3.1 Computational Results 

3.1.1 Graphite 

We started out with a study of bulk and monolayer graphite'. The problem 
with the ASA-approach is that it is only applicable 1 to systems ol reasonably 




close-packed spheres. (Jraphite is a very open structure as is evident Irom 
the weak coupline: between the layers and the rather open honeycomb lattice 
in the 2D-layers. The system can nevertheless be described in terms ol close- 
packed sphere assemblies by introducing a targe number ol empty spheres. 
By placing empty spheres in the hexagon centers ol a graphite layer, one 
produces a closed packed layer ol spheres. ( lose-packed stacking of similar 
I avers, some completely consisting ol empty spheres and some with carbon 
and empt v spheres allows us to build up the graphite st met ure. Bulk graphite 
with the AB stacking is obtained by inserting two layers of empty spheres. 
The graphite layers become almost completely decoupled when three 1 layers 
of empt v sphere's are inserted. This allows us to study both a monolayer and 
the* actual crystal structure of graphite. We can slighlty change' 1 lie 1 size of 
the empty spheres inserted between the graphite layers and thus modi I v eja 
of t he struct ure. 

The band structure's obtained at the experimental lattice constants for 
monolayer and bulk graphite are shown in figs. 1-2. These results are 
in excellent agreement with previous calculations in the literature [I. «>. (>] 
including the weak splittings of the bands in induced by the coupling ol the 
layers in the bulk graphite crystal structure'. 

Fig. :} shows the monolayer density of states (DOS), which is very similar 
to that of bulk graphite. Note the position of the fermi level at the DOS 
minimum in between the x-bonding and anti-bonding peaks, related to the 
critical points al M. The high peak in the range —10 to — * e\ corresponds 
to a region of overlap between the rr and t bands. 

[ be total energy calculations revealed a problem with the ASA. It was 
impossible to obtain a meanignful c/a ratio and even when c/a was fixed to 
the experimental value the theoretical equilibrium lattice constant a was 10 
% larger than the experimental value. Correspondingly, the cohesive energy 
was :}(i Vv too large. Typically for semiconductors or close-packed metals, 
the ASA yields the lattice cosntant to within a few % and the cohesive 
enrgies arc* ~!() { A overestimated. We atribule ihis failure to the inadequate 
description of the charge density tails inside the empty interstitial region. 
The ASA makes a spherical approximation to the latter inside each atomic 
sphere 1 . Since' these (Targe densities decay to almost zero over the length 
scale of a sphere', there must be strong dipoles which have not been taken 
into account. This misrepresentation ol the intersitital charge densitiy in 
turn affects the bonding in the layers. Similar problems are known to occur 
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ior surfaces (see next section on (hi). A possible solution to this problem 
would be to incorporate the dipole moments ol the charge densities inside 
the spheres. The program modifications required to do this could not be 
acomplished within the limited timeframe of this grant. A better solution, 
would be to use the full-potential version of the LMTO method, in which 
charge density and potential are allowed to be completely general. Programs 
to do tliis are already available. 

We decided nevertheless to continue the ASA investigation with (lie more 
limited goal of providing insights in 1 lie electronic structure. File latter 
appeared to he adequately described by the approach. 

3.1.2 Copper 

In preparation to t he calculation of ( 'u/graphite interlaces, we needed to 
calculate Cu ami the (hi surface with a similar stacking of interstitial empty 
spheres as used for the graphite calcualt ion. This facilitates the description 
of the interface system since this allows the same assembly of (dose-packed 
and almost equal size spheres to be used for the various systems. 

The test for bulk (hi total energy properties (equilibrium lattice constant, 
bulk modulus, cohesive enrgy) yielded excellent agreement with a standard 
calculation in which no empty spheres are used. 

Next, we proceeded to the surface calculal ions. A supercell ot (i layers ol 
(hi separated by empty sphere layers was used. II ie resulting surface and 
bulk-like (center of (hi-slab) densities of states (DOS) are shown in Fig. I. 
They show the well known narrowing of the surface FDOS compared to the 
bulk DOS which is due to the reduced coordination. Again, however, the 
surface energy y s defined as the total energy difference between the surface 
slab and an equivalent number of hulk atoms, is very much overestimated. 
These results are similar to what is obtained using a more conventional sphere 
packing. 

From the work of Skriver and Rosengaard [7]. we know that this prob- 
lem is essentially related to the neglect of the intra-sphere dipole moment 
perpendicular to 1 lie* surface. As for graphite, the results are satisfactory for 
a qualitative study of the electronic structure, but not for a quantitative 
evaluation of surface energetics. 




3.2 Cu/graphite interface 

In order lo si inly t ho ( 'u/graphite interface. we employ a supereell consist ing 
of 1 layers of Cn in a licp slacking and 1 layer ol graphite. 1 his stacking is 
repeated periodically and describes the essential local bonding of the 

('infer {111}|| (bgraphite {0001} 

configuration. (The reason for choosing a hep Cu stacking is that we need 
to have a reflection symmetry in the graphite layer in order to have two 
equivalent interfaces and want to reduce the number of atoms per cell. W it h 
just 1 Cu layers, the easiest solution was to stack them in hep order. This 
has no elfed on our conclusions. I 

One laver of empty spheres is placed bet ween Cu and graphite, leading 
to an ideal structure with a Cu C distance 1 of 2.08 A as a starting point 
for our investigations. This distance is ~U.l A smaller than a typical M ( 
distance in transition metal carbide's, it is thus expected to allow ronsider- 
able interaction between Cu and graphite. The C C distance in the plane 
is 1 . fT A. slightly larger than that in graphite 1 (1.12 A). Nevertheless, the 1 
band structure ol this graphite monolayer is still e-'ssentially that e>l graphite, 
characterized by strenig er anel T-bonding. 

The layer projected densities e»f states (DOS) oblaiued in this niex I<*1 are 1 
shown in Fig. ■">. The Cu anel graphite DOS are essentially unpert urbeel 
from I lie 1 corresponding bulk DOS. Furt hermore. the 1 Cu A*/ baud lie's we'll 
separates! Irenn the majority ol C states, billing in the 1 lenv DOS region near 
I lie* Fermi energy. 1 his is clearly inelicat-i\'e ol a we>ak interae't ion be i twe , eu 
the two. Cpon closer inspection, euie may notice an upward shift by -~0.7 
e V of the position of the Cu :h/ band peak with respect to the Fermi leve'l 
on the interface layer. The Fermi level on the graphite monolayer ran be 
seen te> lie ~ I .(i eV higher than in the isolated monolayer, where it occurs 
at the minimum in the DOS indicated by the arrow in Fig. •">. This is 
clearlv indicativi- of a charge transfer from Cu to the graphite. Because of 
the empty spheres, it is not straightforward to assign tin 1 charge transfer 
qiiantitavely on the basis of the calculated sphere charges. From the area 
under the FDOS curve between the Fermi level position in bulk graphite and 
at the ('u/graphite interface, we estimate that ~U.l elect ron/( '-atom have 
beeen transferred to graphite. Notice, however, that in our model there is a 
Cu layer on each side of the graphite monolayer, for a monolayer absorbed 




on ('ii. we woiil<l thus expert a smaller charge transler. We note here that 
for hulk graphite in contact with hulk metal, the hulk Fermi lewis must, of 
course, he aligned. I lie local shilts ol the korun level near the inteita<e seen 
here, must thus he interpreted as an opposite' shiit ol the hands (and thus 
t he local potential) with respect to the fixed Fermi level. 

We note that there exists experimental evidence lor metal to graphite 
charge transler in other metal systems. (See next section.) I ntort unatel\ . 
we could not determine quantitatively how this ailects tin 3 total energy of 
bonding (i.e. the adhesion energy). 

In the calculations described above tin* ( u was placed directly above 
tin 1 (' atoms of the graphite. An alternative position would he to place 
it over the holes in the honeycomb structure. I urt hermore. the distance 3 
between (' and (hi should be adjusted to minimize 1 the energy. Because 
of the inaccuracies encountered in ihe ASA calcualtion ol total energies lor 
graphite and (hi surface's, we did not atlempl to do so. Ihe interatomic 
distances of our idealized sphere 3 packing probably place us in the' ivpulsiw 
range of the interaction. Expanding these distances, howewr. we change the 
overlap between the sphere's too much to allow any meaningful comparison 
between their total energies. We did a calculation for another (hi position, 
and found qualit at ively similar conclusions about weak bonding and charge' 
transfer. The details of the' DOS near the Fermi level we're' found to de'peml 
on the' del ails of the structure. Since we could not opt imize t he 3 st met tire. we 
cannot- attach too much impe>rtance to a detaile'el analysis ol these features 
in t he' elect ronic st ructure. 

3.2.1 Discussion and Conclusions 

The electronic structure results show that the bonding between (hi and 
graphite is weak. Basically, no chemical bonds are' lonned as is indicated 
by the' FDOS figures. I he' re'ason for this is that the' graphite' monolayer is 
almost- like a c|e>sed shell system. The strong T-bonding related to the cle>se 
inte'rael omic distance's betwevn carbon atoms in the laver le'a\’e 3 s no elec- 
trons free in "dangling bonds to couple to a metal. 1 he (hi d-baud itse'lf is 
practicallv filled. Tims only the (hi Is is available' lor making bonds. Ihe 3 
interaction results basically from a weak overlap between ihe* charge 1 deuisity 
tails. Basicallv the repulsive' e-ore-core interactions dominate and only a wry 
we'ak \ an-ele'r- W aals-t ype of bonding results. I here is some' weak coupling 
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to I Ik* 7r-elertron system which results in some charge transfer. Mils may 
have souk 1 effect on the bonding strength. 

3.3 Literature Study 

3.3.1 Graphite surface energy 

Since* graphite is a layered material with weak interlayer bonding, it is obvious 
that the surface energy of 1 he basal plain 1 must lx* very weak. Basically, 
the surface energy y s in this case is half the interlayer coupling in bulk 
graphite, properly normalized per unit area. Several studies were dedicated 
to calculating the interlayer coupling in graphite*. 

The values are summarized in fable 1. 

Fal)le I: Surface* energy of graphite* {DUO I } 


Kef. V 

BrennanfS] 70 

I) i \ incenzo et al. (!)] 330 

Jansen and Freeman [(>] 250 

(mod et al. [I I] 120 


• Brennan [S] using lattice summations of semi-classical i nt eral exilic 
forces based on quantum chemical calculations, obtained a value* olr 
y s of 70 m.)/nr. 

• I he work of DiVincenzo et al. [9] is based on density functional theory. 
It estimates the interlayer charge density by superposition of the charge 
density tails of self-consistent ly calculated monolayers and treats the 
kinetic' energy by means of gradient expansions instead of using a K oh 11 - 
Sham calculation. Iliese authors, in addition showed that the* bonding 
between tdie layers is rather structure* insensitive since* similar values 
were obtained using a 2D averaged model ( barge density. I hey give a 
very interesting account of the relations between the < dassical Van der 
Waals interaction to the density fucntional theory. I he value obtained 
for y s is 330 iiiJ/in*. 
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• 1 die most arcuate calculations of bulk graphite we're performed i >\ 
Jansen and Freeman [()] using LAPW and by d in and ( oben [ID] usinu 
tin* norm-conserving pseudopotential plain' wave method. Jansen and 
Freeman provide an estimated value lor the energy difference between 
t lie monolayer end the bulk crystal, from which we extract a value for 
the surface energy y 5 = 250 mJ/niA This number, however, is based 
on a coniparsion between calculations perlormed with slightly different 
LAPW codes. (The somewhat older film-LAPW version was used for 
tin* monolayer calculation. ) 

• An experimental value for I he surface energy was obtained by ( Jood et 
al. [U]. These results are based on measuring the heat ol immersion 
in various organic fluids whose' surface energies are known. 1 he value' 
obtained was 120 mJ/mA 

l)IS( VSSION: The DFT values seem to be larger than t he experimental 
value'. Density functional theory is known to ove'rest imate* cohesive energies 
and thus is likely to overestimate the surface' energy by ~ 1 () /A There are 
probably rather large uncertainties also on the' experimental value and ihe* 
older semi-empirical calculation of Brennan. \\ e can probably salely say t hat 
y* < 0.2 ±0.1 J/mA 

3.3.2 Copper 

Several ca Iculat ions e'xist of the surface' emerge ol Copper. I he' value's are 
summarized in I able 2. 

Fable 2: Surface energy of Cu {111} at T = 0 


IM. 

* s (m.I/m 2 ) 

Appelbaum Uamanii [12] 

2016 

Smit h A Banerjea [ 1 3] 

1 S.'IO 

Skriver [1 1] 

1 007 

Miedema [15] 

1 S2”) 

Wawra [l(i] 

2000 




• Appl ebaum and liamann [12] used sell-consistent slab calculat ions with 
an atomic orbital basis. 

• Smith and Banerjea[l3] used the so-called equivalent crystal theory. 

• Sk river [1 1] used the LMTO surface (Teen's fucntion method. This 
approach is closely related to our own ASA-LMTO calcualtions. but 
includes the < , ssential surface' dipole contributions. 

• Miedema [ 1 o] used a semi-empirical approach based on relations be- 
t ween vaporization energy (cohesive enrgyl. surlace energy and elec- 
I ron density. The numbers cited by him are bascially obtained from 
liquid drop measurement s ext rapolated to zito temperat ure using scal- 
ing relations derived in his papers. 

• Wawra's [Hi] approach is also semi-empirical and is bast'd on an as- 
sumed proportionally between surface energies and elastic constants. 
Measurements of the elastic cosntants as function of temperature allow 
him to predict the temperature dependence o( the surface energies. 

Wawra [Hi] gives a value of about 1 ‘500 mJ/m J for the liquid metal at the 
melting point and I (>()() m.j/iid for the solid at the melting point. 

( T)N( ’LUSION: There appears to be agreement on a zero temperature 
value of y s ^ 2 ±0.1 .J/m\ and a value of about l.o ±0.2 .l/m J at the 
melt ing point . 

3.3.3 Analysis of wetting studies 

In the wetting studies of De Vincent et al. [I] a contact angle 0 
measured. Using 

— ! ■? r T’u 

with t he interface energy and 

11 — ;('u ± ");;?■ ;<//■/ ( ’ }i ' (— ) 

— v ui 1 ± eos 0 ). (•? ) 

and the values deduced above for the surface enrgies of (hi and graphite, 
we obtain the adhesion energy 11 to be ~0.2 J/m J and the interface' energy 
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snn'IA) .I/nr. ( liven t he uncertain! ies. we (’an essent ially take* t he* sitrfan 1 hi- 
rtvv of graphite and the' adhesion to ( 'u 1o be zero, whirl) makes the interface 
energy about- equal to the surface enrgy of (hi which is about 2 J/mA 

Since the experiment relates to a liquid metal, one should more properly 
subst it ute a value of about 1.3 .J/m J insterad of 2 J/nr in the* above. In their 
analysis. De Vincent e>t al. [1] oblained rat her different values mainly because 
t hev use a higher value lor the suriace energy ol graphite. A higher value of 
the surface energy of graphite may he expected if the physical surface in the 
(experiment, corresponds to an average over the basal plane and other (higher 
energy surfaces of graphite) due e.g. to t he existence of steps on tin 1 suriace. 
Their value (0.8 .J/m 2 ) is based on the work of Mortimer and Nicholas [IT] 
which however, is based on the questionable assumption that the interface 1 
energy between a metal and graphite is the same as lor the metal with the 
corresponding carbide. We believe' this to he an overest imateol the suriace 1 
energy of graphite. 

This doe's not. of course alfect the 1 validity of DeVincent et al.'s measure- 
ments of the contact angles as function of alloy additions. Contact angle' 
measurements on IIOPC graphite may be useful to obtain a bet ter e*st imate 
of t he true' contact angle on the basal plane' and hence* the present estimates 
of 1 h e a< 1 lx *sion ene-rgy. 

3.3.4 Other studies of metal/graphite interaction 

Several studies exist of the interaction between graphite thin films or mono- 
I a \ (ms and metal surface- 1 *. One* context in which these arise* is in studies of 
catalysis. Deposition of carbon containing molecules (CO. Cll r . e*(c.) <m 
metal surfae'es leads to the formation of a graphite monolayer on t he surface 
under suitable heating treatments. Essentially this occurs under conditions 
of a sufficiently high carbon ('overage which allows the carbon atoms to get 
close* enough to e*ach other to form a - bonded system. In the* case* of a lower 
coverage, the carbon atoms prefer to make a stronger bond to the substrate, 
which leads to a suriace phase* resembling the* bonding in carbide's and there- 
fore called carbidic. A good summary of the physics of these surface* systems 
can b<* found in the theoretical work of Feibeiman [18]. 

Several experimental sludies ol this also exist, e.g. work by Hosed et 
al. [lb. 20]. Papagno e*t al.[21], McConville et al. [22] on Ni/graphile and 
Aizawa e*t al .[23] on TaC/graphite. Several of these authors mention 1 he* 




existence of a charge transfer between 1 he metal and graphite as evidenced 
I ) v pholoemission and EELS. The work of Aizawa el al. [22] showed that 
tins can lead to a softening of the bonds in the graphite layer as evidenced 
by phonon softening measured by EELS. 

In t he case of Ni/graphite t he interact ion was found st rong enough t o force 
a registry of the graphite monolayer with the Ni substrate, as evidenced by 
SEX AES [20] measurements. 

We also found evidence for a charge transfer between Pi and graphite in 
I he photoemission studies ol Eppell [21] al ( A\ HI . 

Finally, we mention that graphite metal interactions have been ol interest 
in the context of scanneling tunneling microscopy, (daphitr is a very much 
used substrate in these studies and metal tip to graphite interactions have 
hence been studied. W e mention for example t lie 1 heorel ical work ol 1 omanek 
and Zong [22] using first-principles density tunctional calculations ol a metal 
atom interacting with a graphite surface. 1 h is work shows that the repulsive 
interaction is weaker for an atom approaching the surlace above the hollow 
position than for the metal on top ol carbon position. I lie bonding minimum 
occurs at a distance of ~ 2 . (> A and is extremely weak. This confirms our 
qualitative conclusions ol weak bonding. A semi-empirical paramet riznt ion 
of tin 1 bondign curve 1 obtained from first -principles is given in this work. I his 
may be useful for future 1 simulations of Pd/graphite interaction. 


4 Conclusions and Recommendations 

Our computational results indicate that the bonding between (hi and 
graphite is weak because of the "Hosed shell’* nature of graphite 1 and the 
full el-shell of (hi. We nevertheless found some indication ol charge transfer 
which is worth further study. Also, transition metals with partially filleel 
el-bands may show somewhat stronger interactions with graphite 1 . Ill is is 
indicated by the work on monolayers ol graphite on these metals. I he en- 
couraging improvements in wetting found by alloying additions ol transition 
medals ( I M ) I o ( u also point in t hat direct ion. 

An obvious line* ol future work would thus be to (extend the work to 
transition metals. In order to obtain accurate quantitative results on llie 
e'lKM’uet ics. full-potential calculations are* required. 

An important question is whether a substantially stronger bonding oc- 
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nirs lor t lies e I'M metals at abrupt i nt <*rlar< % s with graphite. or. whether the 
1 1 n p i*( >vet i bodnmg is due to carbide format ion and < heim< al leactioii. 1 he 
wet t ing st udies point t o t he existence of an i lit er facia l phase 1 and chemical in- 
t('i-act ions, it is of great interest scient ifically to st udv t In* early st age's of t his 
process. The transition metal carbides have tin' rocksalt structure. \ iewed 
along the {111} direction, this structure 1 consists ol alternat ingly carbon and 
TM atom layers. 'The only difference 1 with graphite/ I’M superhit t ices of the 1 
tvpe we have studied is the number ol carbon atoms per I M atom in the 
carbon laver. Transition metal carbides have' very strong bonding resulting in 
refractory compounds. There must thus be a gradual transition in bonding 
strength from the metal/graphite interface to t lie carbidic bond, hssential 
lor the formation of 1 lie' carbide appears to be 1 that the carbons remain suffi- 
ciently far apart so as not to form T-bonds among each other. 1 he < p lest ie>n 
is thus how the metal inserts itself in the 1 first layers ofgraplitie and breaks 
up i he 1 strong graphitic interlayer bond. Studios on model systems could 
provide important insights in the way t he carbide 1 interfacial layer is lormoel. 
Obviouslv. diffusion anel kinetic effects will also be of groat importance in 
t his react ion. 

ft veil when we accept the formation of a transition metal carbide. 

1 he question remains how the latter bonds to graphite 1 . At an abrupt 
graphite/ TM-carbidc interface, the bonding might be similarly weak as for 
I In 1 graph! ie/ met a I interface. The imporved bodning may be due to some in- 
terlocking of phases and bonding across lion-basal plane 1 surface's of graphite 1 . 
This aspect could be 1 studied by IMectron-microscopv observat ions of the in- 
terfacial layer. 

We also recommend surface science studies of the interaction between 
graphite and metal by depositing a thin layer of graphite on various 1 ransit ion 
metals (or vice versa.) as a function of temperat ure. We also recommend doing 
the contact angle measurements in THY and on Il()P(! graphite in order to 
avoid any extrinsic effects. 

ftinallv. with continued support for our theory project, reliable values 
for the energetics could be obtained using full-potential calculations. I he 
ftP-LMTO method tins been proven to provide accuracies comparable to 
the LAPYY metliod. The latter has proven 1o provide accnrale values tor 
graphite and metal surfaces. Future theoretical research projects related to 
the present work could include: 




• full-pol rut ial LMTO calculations of graphite, ( ’u-surlace energies and 
( ’u/graphite interfaces: 

• similar studies tor a transition metal such as li and ( r: 

• st udies of systems intermediate between Fi/graphite and I i( . 

Accomplishment o( these research goals however. will requite a moie sub- 
t ant ial elforl than was possible withing the framework ol this grant. 
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ficfrf captions 


Fig. I Hand structure of a. monolayer of graphite. 

Fig. 2 Hand structure of AB-stacked hexagonal graphite 

Fig. i Density of states of a monolayer of graphite 

Fig. 4 Local densities of states of (hi {111} surface 

Fig. o Local densities of states at Cu { 1 1 I } /graphite {0001} interlace 
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